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Detonation in a tube with an open space between charge 


and tube wall 


By Henrix L. SELBERG! 


With 5 figures in the text 


In the present paper the detonation of a cylindrical charge lying in an un- 
deformable cylindrical tube is studied theoretically (Fig. 1). Charge and tube 
are assumed to be co-axial, the charge being initiated at the one end. The charge 
detonates with a velocity which is assumed to be constant. After passing through 
the reaction zone the gas usually flows backwards in relation to the front, but 
in certain cases a flow is observed ahead of the front in the space between 
charge and tube wall. In the latter case a gas front is generated which moves 
with a velocity greater than the detonation velocity of the charge. 


Symbols 


x =distance from end of the tube where the initiation takes place. 
t =time, ¢=0 at the moment of initiation. 
D =detonation velocity of the charge. 
D,=front velocity of the gas. 

Fg =the gas front if Dg>D. 

é =2/ tD: 

§g =D,/ D. 

p =pressure. 

@ =density. 

eel / 0. 

E =internal energy per mass unit. 

S =entropy per mass unit. 

c =V (dp/ d@)aq sound velocity. 

u =particle velocity. 

@ =cross-sectional area of the charge. 

A, =internal cross section area of the tube. 
A,=A,~—a. 


Index 0 indicates the intact explosive, index 1 the gas at a cross section 
L, of the tube at a constant distance /, behind the reaction zone and index 2 
(if D,>D) the gas at a cross section L, of the space between charge and tube 
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wall at a constant distance J, ahead of the reaction zone. In L, and L, velo- 
ties, pressures, etc. are assumed to be constant over the entire cross section. 
Further, we assume that the reaction takes place in the space between L, and 
L, and that dS=0 after the gas has passed these cross sections. D indicates 
the detonation velocity for the explosive at the charge density 0)=0)4/A, ac- 
cording to the Chapman—Jouguet condition. Index g indicates the state at the 
gas front F,, and index B the state at the end of the tube where the initia- 
tion has taken place. 


Fundamental equations 


We shall now study the parts of the tube where the entropy remains con- 
stant. The time ¢ is assumed to be so large that 1,/Dt and 1,/Dt may be 
ignored. In such circumstances the state of the gas must be function of € and 
we obtain for example for the pressure [1] 


dp(dx xdt 
Loaavng 


and from this follows dp=0 if at least one of the conditions 


dp dx 
d é =0 and dt = ED 
is fulfilled. 


Thus the pressure propagates with a velocity &D in every point & for which 
dp/dé+0. 
According to Riemann, on the other hand, the quantities 


Dp D 
d 
praut [2 and pau [2 
oc oc 
will be propagated at the velocities w+c and u—c. 


The detonation gas before the reaction zone 


_Here py is constant. The pressure p is consequently either constant in the 
vicinity of € or y is propagated at the velocity 


u—c=ED. 
284 


ARKIV FOR FysIK. Bd 18 nr 17 


y &3 E4 Eg 
Fig. 2. Fig. 3. 


If p,>~p, the pressure distribution ahead of the front will therefore be the 
following: 
In the interval 1<&<é&,, in which &, is determined by the equation 


the pressure p remains constant = Py. 
In the interval §,<&<€, in which &, is determined by the equation 


Ug — Cg = &, D, 


the formula uw—c=£€D applies. Here 


De 


want [ 2. | (1) 


dD 


In the interval §,<&<6,, finally, p remains constant = p,. 

If p,=p., then p=p,=~p, holds good within the entire interval 1<&<&,. 

Here we have assumed adiabatic expansion. It is clear, however, that if the 
gas flow generates a shock wave front at ¢’ where 1<&’<&,, then this shock 
wave will proceed backwards at supersonic velocity until it disappears. p,> py 
therefore always applies. Fig. 2 shows schematically p versus & ahead of the 
front if p.> py. 

If we assume the interspace between charge and tube wall filled with air 
from the beginning, then p, can be determined with the aid of the Hugoniot 
curve of the air. In Fig. 3 the unbroken curve indicates pressure versus piston 
velocity for shock waves in air. The point of intersection between this curve 
and the curve of dashes representing eq. (1) gives the point pg, wy. 


The detonation gas behind the reaction zone 


Let us now proceed to the section 0<&<1 behind the reaction zone. As we 
have seen above, the pressure will propagate with the velocity §D in every 
point £ for which dp/d&+0. On the other hand, the propagating velocity is 
according to Riemann equal to u+c, where wu is determined by the formula 
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Pp 
dp (2) 
oc 


Py 


U=Uy, + 


Consequently ut+c=ED 


applies in every part of the interval 0<&<1, where the state is not constant. 

Let us now assume that D,=D, so that no gas flows ahead of the reaction 
zone. We assume that the charge detonates at its normal detonation velocity 
according to the Chapman—Jouguet condition. The state at Z, then corresponds 
to the one obtained immediately behind the reaction zone if the same explosive, 
but with the density 6,=0,4/A, detonates at the detonation velocity D. The 
detonation velocity corresponding to the density @) according to Chapman—Jouguet 
is <D, i.e. the detonation velocity D is in reality supernormal [2]. Behind the 
detonation front a section with a constant pressure extends from the front §=1 
to €, determined by the equation 


U,+c,=&, D. 


Behind £, the pressure falls monotonously with decreasing € according to the 
formula 


u+c=éED 


in which uw is determined by (2). This takes place until w or p attains the 
boundary condition which applies at =0. If this occurs in the point &, then 
U=Uz, P=Ppz, etc. applies for O<E<&,. Fig. 4 shows schematically pressure 
versus € when the tube is closed at the end 2=0. 

In the Hugoniot diagram Fig. 5 for the explosive of density 0,=0,a/A, the 
detonation point with pressure p, and spec. volume V, will be represented by 
the point A. The straight line from the initial point to A corresponds to the 
detonation velocity 


1 p 
D=— = ; 
00 1/65— V; 


Conditions for D,>D 


If the charge detonates at its normal velocity according to Chapman—Jouguet, 
D, will be =D, i.e. the detonation gas will not flow ahead of the detonation 
front. This holds good as long as D>D. But if D<D this corresponds in the 
Hugoniot diagram (Fig. 5) to the detonation point lying on the line of dashes 
which does not cut the Hugoniot curve. The only possible explanation is that a 
part of the detonation gas flows forward ahead of the detonation front. 

This phenomenon is observed when dynamite detonates at a low velocity. 
For 25 mm cartridges this is ca. 2.300 m/s, whilst D possibly falls below this 
value only for very small values of a/A,. The detonation process has been studied 
in various tubes at the Nitroglycerin AB, Physical Research Laboratory [3]. 
It has been found that the detonation fades out in certain tubes, probably due 
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The Hugoniot curve 


0 86 Es 7 v— 7-7, 


Fig. 4. Fig. 5. 


to the dynamite being dead-pressed by the compressed air and gas layer ahead 
of the detonation front. 

The question if a stable, stationary detonation process can develop in an 
undeformable tube if D<D is probably merely of theoretical interest. We 
simplify the problem essentially by assuming that the charge lies protected against 
the pressure from the gas and air layer ahead of the front. Let the released 
energy per mass unit of the detonation gas which flows through the cross- 
section LZ, and/or L, be Q, and/or Q, and let n, and n, be the fraction of the 
explosive which flows through these cross sections. The following equations 
then apply 


2 a Pi a x P2 
D=n, ((D um) +P) ns ((u Whar alt 


2 ee 2 — py 
ee (eo th) +AR,+-9,] +g (S257 az, +29). 


The first three equations express the constancy of the mass, the fourth that 
of the impulse and the last that of the energy. To this shall be added two 
equations of state: 


E, =F, (P1; 01): Ey =F, (Pe, 02): 


We thus have in all 7 equations and 10 unknown quantities. Of the missing 
equations, one at least must express a property of the reaction. The simplest 
suggestion is that the gas which flows through the two cross sections L, and L, 
originates from the same reaction zone and expands adiabatically after leaving 
it. This gives 
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Pi 


(uy DP -(D-u*=2 | (2) 
OQ ad 


Dz 


The condition that n, is maximum gives 
D—u,=¢, (3) 


which is identical with the rarefaction setting in just behind the cross section 
L,. The point p,,u, must lie above the unbroken curve in Fig. 3. If the cal- 
culated point ,,u, lies below this curve, it implies that D—u,>c,, and instead 
of (3) we get the condition that the point p,,u, lies on the unbroken curve 
in Fig. 3. 

A numerical calculation of the possible ways of propagation according to the 
principle described would probably be rather comprehensive and for various 
reasons of no great practical interest. 
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